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Solutions  of  the  hyperbolic  equation  of  turbulent  flow  - 
in  the  near  ground  layer  of  air  with  indifferent  strati¬ 
fication  are  investigated.  As  an  example  the  diffusion  of  • 
smoke  at  ambient  temperature  from  a  chimney,  is  examined. 

The  generalised  expression  for  diffusion  is  indicated  for 
cases  with  arbitrary  temperature  stratification  of  the  air.  {  ) 

Since  the  diffusion  of  a  substance  in  a  turbulent 
medium  comes  about  because  of  turbulent  fluctuations  of  wind 
velocity,  the  size  of  which  is  limited  (and  under  the  actual 
conditions  in  the  lowest  layers  of  the  atmosphere  is  rather 
small)  the  propagation  of  the  substance  in  the  atmosphere 
takes  place  with  a  limited  velocity  and  in  a  number  of  cases 
(e.g.,  after  the  release  of  an  instantaneous  point  source) 
the  propagation  front  of  the  substance  can  be  observed. 

For  describing  such  cases  the  common  expression  for  diffu¬ 
sion  having  a  parabolic  character  is  shown  to  be  inappro¬ 
priate.  In  our  earlier  paper  [l]  a  system  is  investigated 
in  which  allowance  is  clearly  made  for  the  final  size  of  the 
turbulent  fluctuations,  and  to  describe  turbulent  diffusion 
we  derive  a  system  of  equations  of  the  hyperbolic  type. 

The  simplest  version  of  such  equations  for  the  case  of 
diffusion  in  the  vertical  direction  z  only,  is  of  the  form 


-  3  a^a  ~  <u)  “  huqt>  ff2,  =  a  fa*  3s)  +  §£uc[2»  *  •  •  (1 ) 


where  qi  and  q2  are  the  probability  densities  for  the  co¬ 
ordinate  z  of  the  diffusing  particles  moving  respectively 
upwards  and  downwards;  t  is  time;  u  is  the  size  of  the  tur¬ 
bulent  fluctuations;  z  the  velocity  components;  a  is  the 
frequency  of  change  of  direction  of  movement  by  the  diffusing 
particles  (frequency  of  "dispersion").  An  analagous  equation 
was  proposed  earlier  in  the  work  of  V.  A.  Fok  B.  I,  Davydov, 
and  E.  S,  Lyapin  quoted  in  reference  [l]  and  also  in  the 
article  by  S.  Goldstein  [2]  published  in  1951,  and  isex- 
tremely  close  in  content  to  Fok's  work  in  1926.  Equations  (l) 
correspond  to  the  description  of  turbulent  movements  in  the 
direction  z  as  the  aggregate  of  "streams"  having  flow  rates 
u  and  -u,  the  diffusing  particle  transferring  from  streams 
of  one  type  to  streams  of  another  type  on  an  average  with 
frequency  a. 

Examining  only  stationary  turbulence  we  may  consider 
the  values  u  and  a  independent  of  time  t.  Introducing  the 
density  of  the  particles  q  *  qi  +  qa  and  the  density  of  the 
turbulent  particle  stream  Q  ■  (qx  -  qa)  we  reduce  equations 
(l)  to  the  form 


is 

at 


15  = 

oz 


o 
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+  2aQ  =  - 


u6zUq‘ 


.... (2) 


According  to  the  first  of  these  equations  it  is  doss- 
ible  to  put  q  =  Q  =  where  ¥(z,t)  will  have  the'  mean¬ 
ing  of  an  integral  function  of  probability  distribution  for 
the  co-ordinate  z  of  the  diffusing  particle.  The  second 
equation  (2)  takes  the  form  of  a  telegraph  equation  (3)ir 


3  6¥  '  6 2y 

uozuaz  "  Ft2- 


■  The  problem  consists  in  the  determination  of  the  co¬ 
efficients  of  this  equation  u (7)  and  a(z)  for  different  con¬ 
ditions  in  the  atmosphere  and  in  solving  equation  (3)  with  ; 
given  initial  and  boundary  conditions. 

Paper  [l]  examined  the  case  of  diffusion  in  the  field 
of  homogeneous  turbulence  (u  and  a  constants)  in  the  lowest- 
layer  of  the  atmosphere  with  indifferent  stratification* 

In  the  present  paper  we  shall  examine  diffusion  in  the 
ground  layer  of  the  atmosphere  with  indifferent  stratifica¬ 
tion  in  more  detail  than  was  done  in  reference  [l];  we 
shall  define  the  determination  of  the  Riemann  function  for 
equation  (3)  in  this  case  and  examine  an  example  of  prac¬ 
tical,  importance,  i.e.  the  diffusion  of  smoke  at  ambient 
temperature  from  a  chimney.  Furthermore,  we  shall  show 
how  to  determine  the  coefficients  u(z)  and  a(z)  for  the 
ground  layer  of  the  atmosphere  with  arbitrary  temperature 
stratification  and  analyse  equation  (3)  for  these  conditions 

1  *  Diffusion  in  a  Thermally  Homogeneous  Ground  Laver 
of  the  Atmosphere  . 

The  stationary  turbulent  process  in  the  ground 
layer  of  the  atmosphere  with  indifferert  stratification,  is 
characterised  by  a  unique  dimensional  parameter  -  the 
friction  velocity  v*,  so  that  every  parameter  of  turbulence 
haying  the  dimensions  of  velocity  must  be  proportional  to 
v*,  the  characteristic  dimension  of  length  must  be  propor¬ 
tional  to  height  z,  and  the  characteristic  dimension  of  tim 

must  be  proportional  to  -^r.  Consequently  in  this  case 

' *  .  i/** 


u  *  \v*,  a  =  A“»  .  ...  (4 


-3- 


where  X  and  A  are  numerical  constants.  The  constant  X 
should  be  determined  from  experiments  as  a  ratio  of  the  mean 
square  value  of  the  vertical  fluctuations  of  velocity  to  v*; 

O  the  experimental  data  in  our  arrangement  shew  that  the  value 
of  X  is  near  to  unity.  The  constant  A  may  be  determined  by 
taking  advantage  of  the  steady-state  solution  of  equations 
•  (2),  which  has  the  form 


Q  -const,  q(z2)  -  q(zi)  =  ^S-logf^. 

X^V*  *  2 


We  note  that  the  coefficient  in  the  logarithmic 

formula  for  q(z)  must,  as  wo  know,  equal  7^  where  x  is  the 


Karman  constant.  Hence  A  =  Introducing  into  equation 

(3)  the  dimensionless  variables 


z  •  •  Xv*t 

h»-  T  “  1 r> 


•  •  •  • 


(5) 


where  h  is  some  characteristic  height, 
equation  in  the  form 


we  may  re-write  this 


a “V  a2?  2e  ay  n 

dZ*  "  dx*  ~  It  |  3t  ~ 


•  where  e  = 


X 

2k 


is  a  numerical  constant. 


We  are  interested  in 


the  solution  of  this  equation  in  the  range  C>0,  satisfying 
.  .  the  determined  boundary  condition  at  the  earth  s  surface 

C  =  0.  However,  it  is  more  convenient  to  consider  not  half 
the  range  C>0  but  the  whole  range  “00<!C<00»  extending  for  this 


purpose  the  field  of  turbulence  into  the  region  C<0  with 
symmetry  maintained  with  respect  to  the  point  C  =  0.  If  we 
X2V* 

assume  that  a  =  2i^^zT,  corresPoncling  coefficient  in 


equation  (6)  is  of  the  form 


Let  us  restrict .ourselves  to  seeking  a  solution  of 
equation  (6),  corresponding  to  an  instantaneous  point  source 
of  diffusing  substance  at  time  t  ■  O  at  the  point  z  =  h. 

The  initial  conditions  for  the  original  equations  (l)  will 


then  take  the  form 


•  qi  (z,  0)  =  ej.6(z  -  h),  q2(z,  O)  =  £2^(2  -  h), 

i 

where  £1,  £2  are  the  probability  of  positive  and  negative1 
directions  of  initial  movement  by  the  diffusing  particle. 

For  example,  particles  of  smoke  being  discharged4  from  a 
chimney  move  upwards  initially  so  that  in  this  case  £1  *  ‘i , 
£2  =  0.  With  the  explosion  of  a  shell  in  the  air,  on  the 
other  hand,  the  particles  as  a  rule  have  no  preferred  direc¬ 
tion  of  initial  movement;  in  this  case  we  should  assume 
that  £i  =  £2  =  The  initial  conditions  formulated  for 

the  function  V  will  thus  be 

?(0  «  V(C,  0)  =  E  (C  -  1),  9(0  =  (§|)T-0  =  (ea  -  ej(c  -  1),'  ....(7) 

where  E(0  is  an  improper  distribution  of  probabilities, 
i.e,,  the  function  is  equal  to  zero  when  C<0  and  equal. to 
unity  when  C^O. 

.Let  us  first  examine  the  simplest  case  of  a  ground 
level  source  (n  =0).  In  this  case  the  function  ?CC,  t)  can 

depend  only  on  the  ratio  of  its  variables  ~  =  Z  =  — 4—,  and 

T 

equation  (6)  takes  the  form 

(t2  -  1)  +  e}||  =  0 

(the  upper  sign  applies  in  the  range  C>0,  the  lower  when  . 
C<0) .  Since  the  substance  diffuses  with  a  final  velocity 
u  =  \v*,  at  the  moment  t  it  will  be  distributed  accordinq  to 
the  limits  -Vv^tSz^Xv^t  so  that  V(z)  can  be  distinguished 
from  a  constant  only  when  Kl<i,  and  it  must  be  that 
¥(-  1  -  0)  =  0.  ?(+  1  +  0)  =  1.  The  first  integral  of  the 
equation  for  ¥(?)  takes  the  form 


.  _ 

d*  r1  ’  52 

cTC  -  i_C2 

r-V 


Cr 

when  S>0, 
when  Z<0, 


o 
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where  cit  ca  are  constants  of  integration  which  must  not  be 
negative  since  the  density  of  the  diffusing  substance 

O'2*1)  =  df  *  is  not  negative .  The  expression  for 

7(5)  whi'-.h  satisi'es  the  v ondii ion  .  when  C  =  +  (1+0)  has 

the  form  “ 


1 

y(c)  =  { 


-  when  0  < ' s  1 » 
when  - 1  s < •<  °. 


At  the  point  5  =  0  the  function  7(0  has  the  discontin¬ 
uity  7(+  0)-7(-0)-1  -  2z~2'  va^-ue  which  has 

the  meaning  of  that  portion  of  the  diffusing  substance 
existing  at  the  moment  t  at  the  earth’s  surface.  This 
value  does  not  vary  with  time  and,  naturally,  may  be  equal 
either  to  unity  or  to  zero.  In  the  first  case  we  get 
Cr  =  ca  =0,  i.e.  7(5)  =  E{z)  which  denotes  the  absence  of 
diffusion.  Rejecting  this  uninteresting  case  we  will 
assume  ci  +  ca  =  2e.  Then  in  the  semi-distribution  £<0  the 

portion  of  the  diffusing  substance  ^  exists  the  whole  time 

Qnot  varying  with  time  and  consequently  equal  to  its  initial 
value  ea.  We  thus  get  ca  =  2eea,  ci  =  2eei.  Finally  we 
have 


7(5)  = 


1  -  ei(y-~|-^J-)  when  0  £  5  £  1, 

62  when  -  1  £  5£  0.  ....(8) 


Here  7(0  is  continuous  at  the  point  5=0  (since 
ei  +  £3  =  l)  but  the  density  of  the  diffusing  substance 
q(z,t)  undergoes  a  perturbation: 


q(  -  0,  t) 


2eeg 

Xv*t’ 


q(  +  0,  t) 


26  6! 
Xv*t  * 


The  existence  of  this  discontinuity  is  reasonable  since 
diffusion  of  the  substance  at  the  point  C  =  0,  at  which  the 

frequency  of  dispersion  a  =  is  infinitely  great,  is 

extremely  difficult.  3y  analogy  the  absence  of  discontin¬ 
uities  in  the  function  ?  (C.-r)  at  the  "fronts"  jC[  =  x  is 
explained,  i.e.  the  equality  to  zero  of  the  probabilities  of 
.finding  the  diffusing  particles  at  points  z  =  +  Xv*t..  In¬ 
deed  each  particle  discharged  from  a  source  2  =  0  almost 
certainly  undergoes  dispersion  while  still  near  the  source 
and  consequently  does  not  reach  the  "front". 


When  ei  =  1,  62  -  0  all  the  diffusing  substance  at  any 
moment  in  time  is  in  the  half  range  z>0  so  that  this  case 
corresponds  to  the  condition  of  "reflection"  of  the  substance 
from  the  earth’s  surface.  Then 


q(z,  t) 


__1 _ 

xv^t 


0  <;  2  £  Xv*t. 


The  character  of  distribution  of  the  substance  in  height 

depends  here  essentially  on  the  number  a  -  ^7,  i.e.  on  the 

ratio  X,  the  sire  of  the  vertical  fluctuations .of  wind  speed 
to  the  friction  velocity.  Theoretically  there  are  two  con¬ 
ceivable  forms  of  "front"  for  the  diffusion  of  the  substance. 
If  X>2x=0.8,  i.e.,  if  s>1,  then  when  approaching  the  "front" 
the  concentration  of  the  substance  tends  towards  zero  (since 

X2  v* 

the  frequency  of  dispersion  a  a  2*  z  ^en  comparatively 

small),  so  that  the  "front"  does  not  appear  as  a  surface  of 
sharp  discontinuity.  But  if  X<2x,  i.e.  e<l ,  then  when 
approaching  the  "front"  the  concentration  increases  (since 
the  frequency  of  dispersion  under  these  conditions  is  com¬ 
paratively  small),  so  that  the  "front"  appears  as  a  surface 
of  sharp  discontinuity.  Under  actual  conditions  in  the 
ground  layer  of  the  atmosphere  it  appears  that  the  first  of 
these  cases  is  realised,  i.e.  the  value  e>1  applies. 


In  the  presence  of  a  continuous  line  source  of  diffus¬ 
ing  substance  perpendicular  to  the  wind  direction,  and  ne¬ 
glecting  horizontal  dispersion  in  the  wind  direction, . the 
formula  for  the  concentration  is  found  by  replacing  time  t 
in  (9)  by  4,  where  x  is  the  co-ordinate  along  the  wind 
direction  and  U  is  the  wind  speed,  considered  constant.  .We 

thus  get 


-7 


q(x,  2) 


(1  -5 

_e_  ). - Ezil —  0^2^  2a£s 

XVWX  (<  +  Ur  \s~ri*  11  u  * 

'  Xv*x' 


The  vertical  profiles  of  concentrations  at  various  dis¬ 
tances  x  from  the  source  are  seen  to  resemble  one  another 
and  their  form  is  characterised  by  the  function 
f  (CO  =  (l  -  £  )e  -1  (*  +  C)-2"*1  the  graph  ox  which  with  values 
e  •=  1,  f-  is  shown  in  Figure  1. 

■  '  2.  Diffusion  from  an  Elevated  Source 

The  solution  of  equation  (6)  under  arbitrary 

inizial  conditions  7( C,0)  =  9  (C)  >(§“r)_=0  =  9(C)  is  supplied 

by  the  Riemann  formula,  which  for  the  given  equation  has  the 
form 


27(C, 


t)  *  <?(C  - 

-eMC, 

•e;  £  - 

0)  +  9(C  +  t)v(C,  %i 

C+-c 

r  yv*  ^ 

r 

T  J 

{v(C 

'  )v  (C, 

•e;  S'.  C)  • 

C-'c 

- 

ax 

(10) 


Here  v(C,  x;  C',  x')  is  the  Riemann  function,  charac-  : 
terisinq  the  influence  of  a  single  “impulse"  applied  at  the 
point  at  moment  x*  which  develops  at  the  point  C  at 
moment  x^'  •  In  accordance  with  equation  (6)  the  distur¬ 
bances  propagate  with  unit  (dimensionless)  speed,  so  that 
the  Riemann  function  differs  from  zero  only  when 
x-x*Z\ C-C'  |  .  As  a  consequence  of  the  independence  of  the 
coefficients  of  the  equation  and  time,  it  (the  Riemann  func¬ 
tion)  depends  only  on  the  difference  of  the  variables  x  and 
x* .  Moreover,  according  to  its  meaning  it  cannot  depend 
on  the  value  of  h,  which  enters  into  the  scales  of  measure¬ 
ment  of  length  and  time,  i.e.,  it  can  only  depend  on  the 
ratios  of  its  variables.  Consequently,  the  Riemann  func-  " 
tion  in  our  case  must  take  the  form 


v(C,  *;  C\  t')  = 


O 
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where  VU,^')  differs  from  zero  only  when  Jc  -  C' |  £  1. 


Moreover,  from  considerations  of  symmetry  it  must  be 
that  V(-  C,  C’)  »  VU,  - 


For  £'  and  x'  the  Riemann  function  satisfies  an  equa¬ 
tion  linked  to  (6)  i.e.  obtainable  from  (6)  by  substituting 
C'f  x‘  for  C,  x  and  changing  the  sign  before  the  last  term. 
For  VU,  K')  this  expression  takes  the  form 


o2v 

oXT~sr 


•  2c  Wrav  ,  r/6v  >1  _  n 
KTA^oC  +  ?  FCv  -  °- 


For  K.  A  K.'  and  x  A  x'  when  C  and  x  the  Riemann  function 
■satisfies  equation  (6)  giving  for  VU,  S')  an  equation  . 
differing  from  that  obtained  above  by  alteration  of  the 
positions  of  S  and  S'.-  The  difference  of  these  two  equa¬ 
tions  has  the  form 


a 2v  o2v 
af*  ~  W 


+  5 


The  required  solution  of  this  equation  is  clearly  determined 
by  the  known  conditions  in  the  characteristics  of  equation 
(6)  i.e.  in  straight  lines  x  -  x'  =  +  (C  -  C')  or  S  -  K*  = 

+  1. 

The  indicated  conditions  for  the  function  K1)  take 

the  form 


dlogV ( K,  L-L-ll  _  _e _ „  dlp--qy js.  =  ,  e_ _ 

dS  -  nr  +  rp  d s  ITTTs  +  ip 

whence  are  obtained  the  following  values  of  the  Riemann 
function  for  the  characteristics 

r(r+r)  when « >  °. 

V(5  +  1,  O  -  v(5.  «  +  1)  4(5_J_1)E  when  <  <  -  1, 

0  when  -  •  £  C  £  0. 

-9- 


....  (12) 


Hence,  in  particular,  it  follows  that  the  function 
V(5,  5  )  is  symmetrical  with  respect  to  its  variables. 

OThe  particular  solution  (3)  found  above  of  ecua .ion  (6) 
immediately  permits  determination  of  the  values  of  tie 
Riemann  function ‘V (5,  S')  when  S'  =  0.  Actually,  substi¬ 
tuting  in  the  Riemann  formula  (10)  the  initial  conditions 
(7)  we  get  after  simple  xe-arrangement  with  5<1 •+ 


7(C,  *}  -  (C  -  %  -  l)v(J,  |-l)+ijE(C+T  -  i)v(J,  J  +  i) 


i  + 


"  ♦  ?)  ♦ »  r  *(<'  -  <0 

T  *  1 

{«v(S.  f)  »v(|,  t')-,  ,  ' 

•  - ST - * - of> — jd<  ♦ 


....(13) 


Proceeding  in  this  formula  to  the  limit  where  h  -*■  0  we 
axe  satisfied  that  V(S,  0)  is  the  item  in  solution  (8) 
corresponding  to  a  ground  instantaneous  point  source  of  the 

admixture,  which  has  the  coefficient  -■&  ^  E-1.  Eliminating 

this  item  in  formula  (8)  we  obtain 


V(c,  0)  =  v(0,  O  =  (j-f-flf)6- 


.... (14) 


The  complete  expression  for  the  Riemann  function  will 
be  sought  in  the  form 


V (5 ,  K')  = 


-  1 


l  UI  +  [5'|  +  i 


W(2),  Z 


•  \z 


=  1,1, (1  -  ill 


455' 


....  (15) 


This  expression  possesses  the  requisite  properties  of 
symmetry  V(5,  5')  =  V(5/,  5)  =  v(-  5,  -5').  By  the  char¬ 
acteristics  5  -  5'  =  +  1  the  first  factor  is  transformed 
into  expression  (12), "while  the  variable  Z  of  the  function 
W  is  zero,  so  that  W(0)  *  1.  The  area,  of  determination 


U,  £  1  of  the  function  V(S,  ^')  naturally  breaks  _nto 
three  parts  in  v.'hich  the  function  W(z)  can  be  determined  by 
difference  (rigure  2.).  In  the  first  place  it  allows  us 
to  distinguish  the  areas  in  which  K  and  K ‘  have  the  same 
signs,  from  areas  in  which  the  signs  of  i  and  x,'  are  diff¬ 
erent.  In  the  second  place,  for  the  characteristic 
K*  -  £  =1  at  points  5  =•  0  and  5  =  -1  the  Riemann  function 
is  not  analytical  and  since  the  properties  are  distributed 
along  the  characteristics,  it  is  possible  to  anticipate 
breaking  of  the  analytical  ability  of  the  Riemann  function  . 
at  tne  characteristics  S '  +  s  =  +  1  which  separate  at  the 
part  of  the  area  where  the  signs ~of  c,  and  I.'  are  alike. 
Therefore  .we  consider  three  zones 

I  1.5  -  5'|  2  1  .S  |5  +  Z'\, 

II  U  +  5'|  £  1,  £  0, 

III  |5  -  5'|  S  1,  |5  +  5'J  £  1,  55'  <:  0. 

Substituting  (15)  in  (ll),  in  zone  I  we  obtain  for  W(z) 
the  equation  * 


Z(Z  -  l)2W'  +  (Z  -  1)(2Z  -  1  )W #  -  eaW  =  0. 


In  addition  it  is  necessary  to  choose  a  solution  of  the 
equation,  satisfying  the  condition  V7  (0)  =1.  In  zone  II 
for  W(z)  such  an  equation  is  obtained  and  it  is  necessary 
to  choose  a  solution  bordering  on  the  first  zone  [Z  ~  l) 
coincident  with  the  solution  for  the  first  zone,  and  more¬ 
over  satisfying  the  condition  \‘i  (m  =  l)  emerging  from  com¬ 
parison  of  formulae  (14)  and  (15).  In  zone  III  for  W(Z) 
the  equation  obtained  is- 

Z2(Z  -  1)W*  +  Z(2Z  -  1)W'  t  e2W  =  0, 


the  solution  of  which  must  satisfy  the  conditions  W(-  ~)  *  1 
and  W(z)  =  0(Z"e)  with  Z  -  0.  Using  the  formulae  for  the 
analytical  continuation  of  the  hypergeometric  functions  we  ■ 
obtain 


! 


IO 


(1  -  z) zr i (c ,  1  +  e;  1;  Z)  in  zone. I, 

rr2d  +  c)  t  r(i  ■  25)1/2  ne  _  /C|  e;  1  +22-  Z_=J.\ 

Lrii  “2T/  2  r^rr  ~£TA“z—;  aFiv  *  z  ) 


+  rrejS'f rHr^'f”)  aFi(“ e*  - e:  1  -  ^ 

Z)aFi(r  e,  -  e;  1  -  2e ; 


Z  -  r 


zona  II, 


in  zone  III. 
_ (16) 


Here  2F1  (c,  P»  Y»  Z)  is  the  symbol  of  the  hyoergeo- 
metric  function.  Formulae  (15)  -  (16)  specify  the  deter¬ 
mination  of  the  Riemann  function  quoted  in  reference  [l], 
Toqether  with  formula  (13)  they  give  a  solution  of  equation 
(6)  for  the  initial  conditions  (7).  According  to  (13)  this 
solution  with  C<1  -  t  is  equal  to  zero;  at  the  point 
C  =  1  -  t ;  (the  lower  diffusion  "front"  of  the  substance)  it 

/  4  4  V 

has  a  discontinuity  of  value  £2Vi^  -  1,  ~J;  in  the  interval 
1  -  t;  £  C  <  1  +  t  it  is  determined  by  the  formula 

v(c,  x)  •  £v(|,  |  +  1)  +  £a-2-^v(|,  !)• 

1/-T 

while  at  the  ooint  Z  =  1  t  x  (the  upper  diffusion  "front"  of 
the  substance)  it  has  a  discontinuity  of  value  etVl—  +  1,  — ) 
and  with  Z  )>  1  +  x  is  equal  to  unity. 

The  latter  is  demonstrated  in  the  following  way.  When 
C  £  1  +  x  formula  (13)  is  superseded  by  the  formula 


i° 


-  Ml- 1  ■ 

K/x+x 

+  i  c/Jt  (rrrd.  «*)  ; 


0  +  Ml- 1 +  0 

c  av(|.  *')  rAk> 

*  a?  5  HP" — Jd<: 


Differentiating  this  equation  for  <r  and  using  condi¬ 
tions  (12)  and  equation  (11),  v;e  are  convinced  that 

~  °*  At  the  S£:Tie  time  °)  =  1  when  C  *  1. 

Consequently  X7(C,  x)  when  C  £  1  +  x  similarly  equals  unity. 

So  far  we  have  not  touched  upon  the  question  of  extreme 
conditions  for  concentration  of  the  diffusing  substance  at 
the  earth's  surface,  V.'hen  x  <  1  this  question  does  not 
arise,  since  the  substance  issuing  from  the  source  has  not' 
yet  reached  the  earth’s  surface.  The  solution  of  (13)  is 
here  the  probability  distribution.  The  presence  of  rapid 
changes  in  the  function  T (C,  tt)  at  the  points  C  =  1  t  t 
implies  that  at  the  "fronts",  which  are  being  propagated 
from  the  source  upwards  and  downwards,  there  exists  a  finite 
quantity  of  substance:  at  the  "front"  which  is  moving  up¬ 
wards  there  is  a  proportion  of  the  substance  £2(1  -  t)"6, 
diminishing  with  time  according  to  a  power  law;  at  thr 
"front"  which  is  moving  downwards,  there  is  a  proportion  of 
the  substance_  e2(  1  -  x)c,  diminishing  to  zero  with  the 
approach  of  the  "front"  towards  the  earth’s  surface  (pro¬ 
gress  of  the  particles  towards  the  earth’s  surface  is  diffi¬ 
cult  since  the  frequency  of  distribution  is  increasing  here 
without  restriction). 

When  x  >  1  the  solution  obtained  does  not  show  the 
probability  of  distribution  and  consequently  does  not  des¬ 
cribe  in  any  way  the  distribution  of  the  substance*  For 
proof  of  this  it  is  sufficient  to  determine  the  value  7(0,  x) 
when  x  >  1  which  is  completely  specified  bv  the  function  ■ 
V(0,  <;')  =  V(5',  0).  Applying  (13)  and  (14)  convinces  us 
that  YlQ,  x)  -*■  co  with  x  •*  ",  which  proves  the  statement  put 
forward.  It  is  possible  to  shew  that  equation  (6)  has  no 
solution  describing  the  distribution  of  the  substance  for  an 
instantaneous  point  source,  for  which  the  diffusion  of  the 
substance  would  be  through  the  particular  point  K.  ~  0  or  for 
material  accumulating  at  that  point.  This  situation  is 
completely  analagous  to  the  position  for  a  parabolic  equa¬ 
tion  of  diffusion  with  diffusion  coefficient  k(z)  =  xv*z 
which  also  has  no  solutions  describing  the  distribution  of 
the  substance  under  the  conditions  of  complete  or  partial 
absorption  of  substance  at  the  point  z  =  0.  At  the  same 
time  we  can  obtain  a. physically  understandable  solution  of 
equation  (6)  when  x,  >  1  by  introducing  the  condition  of 
reflection  of  the  diffusing  particles  from  the  earth’s  sur¬ 
face.  For  that  it  is  necessary  to  add  to  the  solution  of 
(13)  a  solution  corresponding  to  the  virtual  source  of 
diffusing  substance  at  the  point  z  *  -  h,  i.e,  satisfying 
the  initial  conditions 


qi(z,  0)  *  £2*>(z  +  h),  q2  (z ,  0)  =  ei&(z  +  h) , 
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I 


C 


which  can  be  re-written  in  the  form 

<p(0  =  H{C  +  1),  *(C)  =  (£l  -  £2)6(C  + 


0. 


o  see 
stri- 


Applying  considerations  of  symmetry,  it  is  easy  t 
chat  under  conditions  of  reflection  the  probability  distri¬ 
bution  differing  from  zero  only  when  C  >  0,  should  be  duter- 
mined  by  the  formula  V(C,  t)  -  V (-  *?)  where  '?(£,  is 

the  function  determined  above.  With  x  d  1  this  distribu- ' 
tion  does  not  differ  from  V(C,  x)  but  with  t  >  1  a  diff e- 
rence  will  only  occur  in  the  region  0  £  C  £  r  -  :  in  which 
influence  of  a ^ "ref lected  wave"  is  felt.  At  the  same 

«im+Kthe  01  t«6  ^iffusing  substance  at  the  "front 

the  ref  lec  ced  wave”  ecjuals  zero. 

-  ,-i?he  fesi\lts  obta’f«®d  permit  us  to  make  -calculations 
relating  to  the  solution  of  concrete  problems.  As  an 
example  let  us  calculate. the  concentration  of  the  substance 
u  •  u?  .s  Svif*aCe  with  an  instantaneous  point  source  at 

discharging  all  the  substance  upwards.  Assuming 
that  the  substance  is  not  absorbed  by  the  earth’s  surface 
and  applying  the  preceding  reasoning,  we  must  supoose  (when 
x  2  ij  that:-  ■ 


O 


<1(0,  t)  -  i 


where  Y(C,  x)  is  determined  b/  formula  (13).  The  last 
expression  can  be  put  into  the  form 

•  q(0,  t)  =  5U|> 


where 


t  (?')  =  SXltOju-S.lI  +  W±-^XLl  .  4,r4l  r  5'|)  ’ 

of  o<  -  tr+iTT/s^T" 

After  simple  calculations  we  obtain 


■** 

i 


q(0,  t)  -  %  7^v-rkCr 
h  It  + 


'Til. 
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....  (17) 


V.'ith  h  0  the  seme  result  is  obtained  from  this  form¬ 
ula  as  by  formula  (9),  as  it  should  be. 

Let  us  consider  a  chimney  of  height  h  uniformly  dis¬ 
charging  smoke  at  neutral  temperature  with  velocity  c  g  per 
unit  time.  Assuming  the  smoke  is  finely  divided  (so  that 
its  gravitational  settling  can  be  ignored),  neglecting 
changes  in  wind  speed  U  with  height  and  turbulent  mixing 
along  the  wind  direction  we  may  calculate  the  concentration 
according  to  the  formula 


q(x,  y,  z)  =  cj  6  (x  -  Ut)- - X  ? "tat. 

o  2v'mk~  A  ^  J  • 


Here  x,  y  are  the  horizontal  „o -ordinates  (x  is  direc¬ 
tion  of  wind),  and  k  the  coefficient  of  horizontal  mixing 
in  the  direction  y.  In  particular,  applying  (17),  v;e 
obtain  a  formula  for- the  concentration  from  the  chimney  at 
the  earth’s  surface  in  the  direction  of  the  wind: 


o  nl  -  / >vV*  J_  (s  -  1) r  _  Xv*x  ^  . 

Uh  J  ickh  ir  CC~+T}s'r*»  ?  "  Uh~  *  1  * 

....(18) 

The  abscissa  of  the  maximum  of  this  function  is  pro¬ 
portional  to  h,  and  the  value  of  the  maximum  is  inversely 
proportional  to  h“£2.  Figure  3  shows  the  graph  of  the  uni¬ 
versal  function  -  l)fcU  +  1)“e-1  when  e  =  1, 

characterising  the  form  of  the  concentration  distribution 
at  the  earth’s  surface  in  the  wind  direction. 

Let  us  now  briefly  consider  the  case  where  the  diffus¬ 
ing  substance  is  absorbed  by  the  earth’s  surface.  In  this 
case  the  extreme  condition  c  =  0  should  be  formulated  not 
with  z  0  but  with  z  =  zQ  (roughness  length).  Reckoning 

the  height  z  from  the  plane  zQ  we  must  suppose 
a  =  X-  — — vr-p-p,  so  that  the  function  7 (C,  t)  must  ba  deter¬ 
mined  f r8m  equation: 


a2?  627  2s  '-dV  _  n  r  _  zo  .  n 

oC2,  ‘  6^  ‘ .<£“+W  ot  -  °'  Co  "  Tf  >  °* 


....  (19) 


form 


The  Riemann  function  of  this  equation  will  take  the 


! 


The  determination  of  this  function  and  solution  of 
equation  (19)  unoex  the  given  initial  and  extreme  conditions 
can  be  achieved  in  an  analogous  manner  to  that  stated  above. 

3.  Diffusion  in  a  Thermally  Heterogeneous  Ground 
.  Layer  of  Air 

A  turbulent  system  in  a  thermally  heterogeneous 
_ ground  layer  of  air  is  characterised  b/  three  dimensional 

parameters:  v*,  r^r  and  =r“»  where  q  is  the  vertical  tur- 

c  ^  P  1  ' 

p  o 

bulent  stream  of  heat;  c  and  p,  the  specific  heat  and 

p  • 

density  of  air;  g,  the  accelerating  force  of  gravity;  and 
Tq,  the  normal  temperature  of  the  ground  layer  of  air  (see 

for  example  [3]).  From  the  indicated  parameters  it  is 
possible  to  form  a  length  scale 


In  accordance  with  considerations  of  the  dimensional 
parameter  of  turbulence  u  and  a,  occurring  in  the  equation 
of  diffusion  (3),  it  is  possible  to  write  them  in  the  form 


u 


-**■*(!).■  *  - 


.... (20) 


■•*  where  91(C)  and  92(C)  are  certain  universal  functions.  In 
the  particular  case  of  indifferent  stratification  of  the  air 
scale  L  becomes  infinitely  large  while  the  variable  \  of  the 

A* 

functions  01  and  92  becomes  zero,  so  that  fox  agreement  with 
formula  (4)  (corrected  fox  a  system'  with  indifferent  strati¬ 
fication)  we  must  assume  91(0)  =  92(0)  =  1. 


fr 

f- 

l 

b 


It  is  possible  to  determine  the  function  91(C)  by 
making  use  of  the  equilibrium  equation  of  turbulent  energy 
(see  for  example  [4]J  which  for  the  coefficient  of  turbu-.. 
lence  k  leads  to  the  expression 


O 
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k 


4 


«P. 


where  1  is  the 
Prandtl  sense), 

critical  value. 

square  velocity 

Hence,  applying 


scale  of  turbulence  (nixing  length  in  the 
Hi  is  the  Richardson  nurr.bor  and  Ri  its 

The  value  y  is  proportional  to  the  mean 

of  the  turbulent  fluctuations,  i.e.  **» 

the  condition  91 (0)  =  1 ,  we  obtain 


91(0  «-[l 


L. 


In  order  to  determine  the  function  02(C),  we  examine 
the  steady-state  solution  of  equations  (2),  which  arises 
from  the  conditions 


Q  -  const,  2aQ  = 


-  u67uc*. 


The  right-hand  side  of  the  last  equation  may  be  re¬ 
written  in  the  approximate  form  -  so  that  the  function 

n(z)  changes  significantly  more  slowly  than  q(z).  The  idea 
of  this  simplification  is  to  ignore  the  vertical  stream  of 
diffusing  substance  which  occurrs  when  the  concentration 
gradient  of  the  substance  is  zero  and  which  is  a  consequence 
of  a  difference  in  the  intensity  of  turbulent  fluctuations 
of  velocity  in  different  layars  of  air  (the  presence  of  such 
a  stream  of  substance  is  analogous  to  the  effect  of  thermal 
diffusion).  With  the  indicated  simplification,  by  applying 
(20)  we  can  write  a  steady-state  solution  of  equations  (2) 
in  the  form 


3p  _  0  oc(z/L) 

- - -  y. 


oz 


XV*Z  92  (z/l) 


According  to  the  theory  of  similarity  for  a  stationary 
turbulent  system  in  a  thermally  heterogeneous  ground  layer 
of  air  (see  for  example  reference  [3])  it  must  be  that 


-  St**)  -  -&!>(£■)  -*(?)]. 


where  f(0  is  a  universal  function  which  is  linked  with  the 
Richardson  number  by  the  relation  . 


Ri _ 1 _ 

ii  f'(LY 
*?  1 


Consequently  we  obtain 

: 91(0  *  [1  -  92(c)  =  cf'to  [1  - r^criT2’ 

.-...(23) 

and  equation  (3),  by  calculations  using  equations  (20)  and 
(23),  becomes 


Q-b  -  rrWl  oli1  •  r\zj J  Id  •  If*  - 2£f ' iVll'  -  ridiJ  If  =  0 . (24) 


s.  ~  =  }£2it 


when  the  dimensionless  variables  K  =  — ,  x  = 

4m* 

duced. 


are  intro- 


For  the  analyses  or  thas  equataon  it  as  essential  to 
know  the  nature  of  the  function  f'(C),  We  use  the  equations 
of  motion  and  flow  of  heat  for  a  steady-state  system  in  the 
ground  layer  of  air,  which  have  the  form 

kg  -  (v*)»,  «kg  -  -  -a-. 


Here  U  and  T  are  the  average  velocity  and  temperature, 
and  «<  =  g£ — .  From  these  equations  and  determination  of  the 

CX 

Richardson  nunber  we  obtain  k  =  xv’^L/z /  (C  ).  Substituting 
this  formula  in  (21)  and  determining  1  (in  the  case  of  stable 
stratification  to  which  our  present  discussion  is  confined) 
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we  find 


cm  the  formula 


.  ( .  EL 

1  “  'izl\ 1  ~  Ri  )  » 

/.p 


[as  used  in  the  paper  [4]  (in  which  for  concurrence  with 
experimental  data "we  should  put  p  =  ->)  ]  for  f'(C)  the  alge¬ 
braic  exoression  **  ■. 


r-T^r 


r  •  -  -s/. 

=  C>  i  -  !  /o, 

u  r  _i . 


From  this  equation  it  follows  that  with  C  <  1  then 
* 

^  ^  a  -r  0(C)  and  with  C  “  asymtotically 

f'(C)  ~  1  +  C_5-/3  to  which  the  following  limiting  formulae 
correspond  for  the  functions  (23) 

when  £  <i  1  9i(C)  *=*  1  -  p-C ;  92(d)  ~  1  t  tUC. 
when  C  -*■  03  91(C)  ~  C£/,i2,  92(d)  ~  C*/aB 


Thus  for  small  values  of  C  the  functions  o;  and  92  are 
extremely  close  to  unity  and  the  equation  for  "S  takes  the 
same  form  as  in  the  case  of  indifferent  stratification. 

With  C  —  03  the  fluctuation  velocity  (and  consequently  the 
velocity  of  dispersion  of  the  substance)  decreases  C“3/12 

while  the  frequency  of  dispersion  decreases  as  C  6;  the 
‘  .2 

"coefficient  of  diffusion"  ~  tends  towards  a  constant  value 

2a 

In  order  to  establish  the  parameters  of  equation  (24), 
we  introduce  instead  of  C  a  new  independent  variable 


r.  •  it1  -  rrnJ  dC  *  tv  ~  t1  -  nxil  )  *  t?;1  *  L1  ■  froJ  > 


which  for  small  values  of  C  has  the  fora,  r  ^  C,  and  for 

large  va_ues  or  C  grows  as  C  —  1^.  T we n  ecu 
written  in  the  f or. a 


ation  (24}  v.’ili  be 


^  **» 
c  w  '/ 


o2g 


o  r“  ox 


'1  -  occ(r0—  =  0. 
-  —  S  ^  ox 


....  (26} 


I.'.oxeover,  the  function  g(r)  is  parametrically  repre¬ 
sented  by  the  formulae 


•T: 


g(rt)  =  y-g— ,  n  =  g(l  -  x3^°)  -  ~f(:  -  x--7/^),  0  £  x  £  1. 


For  str.a^l  values  of  r,  g^r,)  =  and  for  large  values 

/ 1  7  \  _L<1T1 

of  r,  asymptotically  g  (q;  ~  A  SraP*t  of  the  junc¬ 

tion  g  (q)  is  shown  in  Figure  4. 

Integration  of  equation  (26)  is  possible  by  the  use  of 
methods  of  approximation. 
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